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Abstract 

The growth and composition of Earth is a direct consequence of planet formation throughout the Solar System. We 
discuss the known history of the Solar System, the proposed stages of growth and how the early stages of planet 
formation may be dominated by pebble growth processes. Pebbles are small bodies whose strong interactions with 
the nebula gas lead to remarkable new accretion mechanisms for the formation of planetesimals and the growth of 
planetary embryos. 

Many of the popular models for the later stages of planet formation are presented. The classical models with 
the giant planets on fixed orbits are not consistent with the known history of the Solar System, fail to create a high 
Earth/Mars mass ratio, and, in many cases, are also internally inconsistent. The successful Grand Tack model creates 
a small Mars, a wet Earth, a realistic asteroid belt and the mass-orbit structure of the terrestrial planets. 

In the Grand Tack scenario, growth curves for Earth most closely match a Weibull model. The feeding zones, 
which determine the compositions of Earth and Venus follow a particular pattern determined by Jupiter, while the 
feeding zones of Mars and Theia, the last giant impactor on Earth, appear to randomly sample the terrestrial disk. The 
late accreted mass samples the disk nearly evenly. 

1 Introduction 


The formation and history of Earth cannot be studied in isolation. Its existence, composition and dynamics are a direct 
consequence of events that occurred throughout the Solar System. Planet formation is a combination of local and 
global processes at every stage of growth. Therefore, the early Earth must be placed in the context of the history of 
the Solar System. 


1.1 History of the Solar System 

There are very few firmly established temporal events during planet formation, however one that seems irrefutable 
is that the giant planets must have accumulated all of their gas envelopes, nearly all of their mass, by the time the 
nebular gas is removed from the disk. Observations of nebular disks about other stars indicate that disk survival times 
are typically only millions of years \Haisch et al.\ |2001[ \Briceno et al. 200 1[ Mamajel^ |2009| l, although there is 
some evidence that disks can last for tens of millions of years {Pfalzner et al. 20141. While we do not have a direct 


observation of our own disk survival time, we do know that chondrule formation-a likely nebular process [Alexander 


et al. 


2008| l ended after a few million years [Kita et al. 2005 [ Villeneuve et al. 2009[ ). Eurthermore, there is evidence 
from modeling of the formation of lapetus including ^^Al decay that it mush have formed ^3.4 to 5.4 My after 
CAIs ^Castillo-Rogez et al.||2007||2009| l. Since the formation of lapetus requires Saturn to be near completion, this 
also constrains the age of Saturn. 

The last giant impact on Earth, the Moon-forming impact, has a date established to be approximately 70-110 
My after CAIs from three lines of evidence: radiometric ( Touboul et al. 2007 Allegre et al. 2008| Halliday 2008Jl, 
dynamical-elemental (Jacobson et al. 2014) l and dynamical-isotopical (Bottke et al. 2014| l. There is some disagree¬ 
ment, others have suggested earlier dates using different models to interpret radiometric data \Yin et al^\20Q2\\jacob-\ 
\sen\\2005\\Taylor et a/.||2009|l, but even these 20-60 My estimates are significantly after the few million year lifetime 
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of the nebular gas disk. Interestingly, the age of Mars is considerably different than that of Earth. Radiometric evi¬ 


dence from the Hf-W system indicates that it likely formed within 10 My {Nimmo and Kleine 20071 and perhaps it 
was within a few percent of its final mass by the time the disk dissipated \Dauphas and Pourmand 201 l| l, so Mars may 
have nearly completed its formation in the presence of the nebular gas. The formation dates of Mercury and Venus are 
unknown. 

The other temporal event with significant evidence in the history of the Solar System is the late heavy bombard¬ 
ment. This intense period of impacts ~400 My after CAIs left a record on the Moon \Tera et oL 1974[ Ryder 1990 
\Cohen et a7[|2000[ Ryder\\lQQ2) , on Earth ^archi et q 7||2OT4| ), and on Vesta through analysis of the HED mete- 
2013 I. Erom the crater record on the Moon, it appears that ^10“^ M^ of material was delivered 


orites {Marchi et al. 

to Earth during this bombardment { Morbidelli et al.\ 2012a[ Marchi et al. 2014)1. Eor a contrasting interpretation of 


the lunar cratering record and asteroid melt chronology, see Chapman et al. ( 2007| l, although this review occurred 
before the terrestrial and HED evidence was presented. 

Nice model There is significant evidence in the orbital structure of the Kuiper Belt that the giant planets migrated 
due to interactions with this disk of outer planetesimals <\Fernandez and |1984) \Malhotra\ |1995[ 


s 112003) 

Levison and Morbidelli 2003 |l. This interaction caused the outward migration of Saturn, Uranus and Nepture and 


the inward migration of Jupiter because the planetesimals are first inwardly scattered by the outer three giant planets 
and then scattered out of the Solar System by Jupiter. When this migration occurred is not directly known, but it has 
been associated with the origin of the late heavy bombardment through a dynamical instability ^Strom et al.\ 2005 | l; 
the ‘Nice’ model stipulates that an early compact configuration of the giant planets undergoes an instability ^Gomes\ 


et al. I 20051. This giant planet instability explains a number of features of the Solar System including the Trojans of 


Jupiter { Morbidelli et al.\ 20051 and Neptune {Tsiganis et al. 2005)1, the number and inclination distribution of the 
irregular satellites of the giant planets ( Nesvorny et al. 2007 Bottke et al. 20101, and dynamical features of the Kuiper 
Belt ^evison et al.\ |2008) \Morbidelli et al. 2014[ ). Eurthermore, close encounters between an ice giant and Jupiter 
reproduce the structure of the main asteroid belt { Morbidelli et a/.[[2010) l. More recently, the model has been updated, 
Nice 2.0, to better reflect the likely orbital spacing of the giant planets in mean motion resonances after removal of the 


gas {Levison et al. \ 20111. 


While a complete rejection of the Nice model would require a number of less likely or so-far incomplete hypothe¬ 
ses to take its place, the connection between the late heavy bombardment and the giant planet instability of the Nice 
model is not required. Eor instance, a Nice model-like instability may occur immediately after nebula gas removal, 
thereby creating the Solar System architecture and small body populations that are observed today. In this alternative 
scenario, the late heavy bombardment would either need to be triggered by the instability of a small fifth planet at ~2 
AU ( Chambers J2007 1, to not occur due to a different interpretation of the absolute age evidence on the Moon ( Chap¬ 


man et al. 20071, or to be triggered by the catastrophic bre ak-up of a Vesta-sized Mars -crossing asteroid due to a 


different interpretation of the cratering evidence on the Moon ( Cuk et al. 


2010 


Cuk 


2012 I. The latter two alternatives 


are unlikely because of new evidence that the late heavy bombardment is recorded in the HED meteorites ^archi\ 

[eF^[20T3) l. 

It is important to note that the Nice model and the Grand Tack model are not the same and are independent of one 
another. The inward-then-outward migration of Jupiter and Saturn called the ‘Grand Tack’ ( Walsh et al. \ 201 1) , which 
is described in Section 2.2 is proposed to take place in the presence of the nebular gas, so much earlier in the history 
of the Solar System than the Nice model. The Grand Tack scenario is consistent with alternatives to the Nice model 
such as the conjecture that the giant planet instability occurs shortly after gas removal and the evidence for the late 


heavy bombardment has some other explanation such as the Planet V hypothesis ( Chambers 20071. 

1.2 Stages of planet formation 


Eor convenience, we divide the history of planet formation into four stages; dust sedimentation and growth, planetes- 
imal growth, planetary embryo growth, and planet growth. While it may be unlikely that these stages are coincident 
locally, it’s very likely that they are occuri'ing simultaneously at different radii across the disk since the orbital period, 
solid surface density, and gas density, viscosity and temperature change dramatically with semi-major axis as well as 
time. Particularly, there is a divide in the Solar System between the giant planets, Jupiter and Saturn, which grew large 
enough to accrete most of their mass from the solar nebula, doing so during the first few million years of Solar System 
history before the nebular gas dissipated, and the terrestrial planets, which did not grow large enough to accrete a 
significant amount of solar nebula gas. Of course, the sub-giants Uranus and Neptune grew at a rate somewhere in 
between, but exterior to Jupiter and Saturn. 

It is computationally difficult (and impossible as of now) for a single numerical simulation to cover all stages of 
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growth, so typically each stage is handled separately. Earlier processes such as dust growth and planetesimal formation 
are often studied in local patches of the disk or in small annuli. This is thought to be acceptable since perturbations 
from elsewhere in the disk are homogenous across the studied zone. Even global simulations of the later stages are 
computationally expensive. 

Each growth stage is marked by a different dominant accretion mechanism for the largest bodies—those we typ¬ 
ically care the most about, but smaller bodies may continue to accrete according to earlier paradigms. We review 
each growth stage, in turn, focusing on new developments, but we begin by defining some new terminology that has 
appeared in the planet formation literature. 

The motion of solid objects in a gas disk is dictated by the Stokes number. It is a dimensionless number with a 
long history in fluid mechanics that quantifies the momentum coupling between an object of approximate radius R and 
the surrounding fluid. In the context of a gas disk, it is defined; St « {ps/po)Ri^/co where ps/po is the ratio of the 
solid to mid-plane gas density, il is the Keplerian orbital frequency, and cq is the mid-plane sound speed. This formula 
is correct for Epstein drag regime, which assumes that the gas mean free path is longer than the object (for a review 
of other drag regimes and the appropriate formulae for determining the Stokes number, see Youdin 2010| l. Accretion 
depends both on the Stokes number of the target and projectile (for recent examples see Chambers 2014[ \Guillot\ 


et al.\ 2014[). When St 1, the object is completely entrained in the gas and when St ^ 1, the object’s motion 


is independent of the presence of the gas. Thus different regimes are identified: objects with Stokes numbers much 
less than 0.01 are ‘dust’, those with Stokes numbers between 0.01 and 1 are ‘pebbles’, and those with Stokes numbers 
just greater than 1 are ‘boulders’. Planetesimals and embryos have Stokes numbers much greater than 1. Erom the 
definition of the Stokes number, an object’s membership in each regime is dependent on the disk properties. The 
Stokes number is the fundamental quantity determining the significance of passive concentration processes, streaming 


instabilities, and pebble accretion; all discussed in detail in Section 1.3 


Dust sedimentation and growth Planet formation begins in the nebular disk that forms simultaneously with the 
protostar from a collapsing molecular cloud as a result of the conservation of angular momentum. Pre-existing dust 
grains as well as condensates from the gas sediment towards the mid-plane of the protoplanetary disk ^Weidenschilling\ 
1980|). Turbulent diffusion of solid material due to entrainment with the gas eventually balances gravity, and so a 


vertical equilibrium structure is assembled in the disk ( Weidenschilling and Cuzzi |1993[|CMZzi et al.\\1993\\Dubrulle 
et al. 1995[ Carballido et al. 20061. In this structure, dust grains can grow via binary collisions accreting through 


adhesive forces, ‘fluffy’ aggregation and compaction, and by mass transfer during fragmenting collisions (for reviews, 
see |Z)om/«/k et al.\\lQYfl\\Iohansen et a/.||20l4| l. Prom these growth mechanisms, a population of pebbles and boulders 
emerges from the coagulating dust. 

Planetesimal growth All objects experience some orbital drag because the nebular gas is partially pressure sup¬ 
ported; a parcel of gas orbits the Sun at sub-Keplerian speeds but boulders orbit at Keplerian speeds. Thus boulders 
experience a headwind, which generates energy loss and their orbits spiral into the Sun ^Whipple\\\912\ . This effect 
is often referred to as the ‘radial drift’ or ‘meter’ barrier, when considering specific disk conditions. Purthermore, as 
objects grow aerodynamic drag decreases, impact velocities increase, and binary collisions are no longer an efficient 


growth mechanism because colliding boulders bounce or even fragment (for review, see Blum and Wurm 20081. This 
is often referred to as the ‘bouncing’ barrier {Zsom et al. 2010|l. 


These boulder barriers restrict growth to S'f ~ 1 or about ~ 0.1-1 m at an AU for a standard model disk \Weiden- 


sizes if ice enhances stickiness and collision velocities remain below ^50 m s 


Wada et al. 

2009 

2013 

Okuzumi 


et al. 2012|l; whether these planetesimals would be too water rich to contribute significantly to terrestrial planet for¬ 


mation is an open question. In Section [T3] we discuss new ideas regarding the interactions between gas and ‘pebbles’ 
that would create ^ 100 km planetesimals directly from <1 m ‘pebbles’ skipping over these barriers Johans en et al.\ 
[20071 [a;iiriF^[2008i i. 

Planetary embryo growth Objects that grow beyond the boulder barriers are planetesimals, and they grow at dif¬ 
ferent rates depending on their size. Smaller planetesimals in the ‘orderly’ growth regime accrete according to their 
physical cross-sections alone and so growth proceeds slowly ^Safronov 1972\\Nalcagawa et al. 19831. It’s possible that 
most large planetesimals including those that continued to grow into planetary embryos and planets bypassed much of 
this growth regime and grew suddenly through gravo-turbulent growth mechanisms from pebble and boulder sizes to 


hundreds or thousands of kilometers across (a description of this growth and its consequences is in Section 1.3 i. 


If planetesimals have grown large enough that the escape velocities from their surfaces match or exceed their 
relative velocities, then gravitational focusing is important ( Safronov and Zy/agi«a[|1969) Greenberg et al. 197811. The 
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enhancement of the collision cross-section due to gravitational focusing is proportional to the square of the ratio of 
escape to relative velocities (for a derivation and further review; Annitage] 2014| l. Since dynamical friction reduces the 
relative velocities of the largest bodies compared to the rest and those bodies already have the highest escape velocities, 
planetesimal growth transforms into a ‘runaway’ growth process; the largest bodies have the largest relative binary 
collision cross-sections and hence grow the fastest and become planetary embryos ( Wetherill and Stewart\ 1989 1993[ 
1992 Kokubo and Ida\\\99b\. While gas drag is less efficient for larger bodies, dynamical friction 


Ida and Makino 


efficiently transfers energy from these larger bodies to smaller bodies that interact with the gas more strongly ( Wetherill 


\and Stewart\\\99'i\^okubo and /t/a||1996| l. Dynamical friction is the damping of the eccentricities and inclinations 
of larger bodies due to gravitational interactions with a large number of smaller bodies. Smaller bodies are also 
more likely to fragment rather than grow due to these higher relative velocities and their lower gravitational binding 
energies, further increasing the number of even smaller bodies and the effectiveness of dynamical friction and gas 
drag ( Wetherill and Stewart 1989| l. 

Runaway growth produces a bi-modal mass distribution in the disk of planetary embryos and planetesimals. The 
formation of planetary embryos is self-limiting and slows down runaway growth by viscously stirring the nearby 
planetesimal population increasing relative velocities ^issauef\ |1987[ \Ida and Making] |1993| l, so growth is slower 
and similar between neighboring embryos. This ‘oligarchic’ growth maintains the bi-modal mass distribution of 
planetary embryos so individually each embryo is about 100 times as massive as the average planetesimal. Orbital 
repulsion keeps planetary embryos about 10 mutual Hill radii apart from each other ( Kokubo and Ida\ 1995 1998| l. 
See \Morbidelli et ar| ( |2012b) l and \Raymond et Qr| ( |2013| l for more mathematical reviews of the runaway and oligarchic 
dynamical processes. 

Planet growth, the giant impact phase Oligarchic growth ends when energy transfer via dynamical friction can 
no longer continue circularizing the orbits of the planetary embryos. Typically, this occurs because the planetesimal 
population becomes too depleted as planetesimals are accreted onto planets and the Sun or scattered from the Solar 
System. Many simulations do not include the role of collisional grinding but this also leads to a reduction of the 
planetesimal population ^Levison et al.| [ 20T2) l. If gas drag was a major source of dissipation, then the sudden removal 
of the solar nebula can trigger the end of the oligarchic regime {Iwasaki et al. 2002[ Zhou et al. |2007| l. 

Without dynamical friction, planetary embryos perturb each other onto crossing orbits leading to either giant 
impacts or scattering events ( Wetherii}\ 19851. The planetesimal population steadily declines during this phase until 
only a few planetesimals are left on semi-stable orbits in the Main Belt. These giant impacts are the last phase of 
planet formation. Once the remaining planetary embryos have found stable orbits and hnished accreting, we call 
them planets. An extended giant impact phase is consistent with the late formation of Earth (^70-110 My; Touboul 


et al. 


2007[ Allegre et al. 2008[ Halliday\ 2()08[ Jacobson et al. 20\A\\Bottke et al. 2014| l. The early formation of 


Mars (^1-3 My; Nimmo and Kleine 2007 \Dauphas and Pourmand 20111 may imply that it hnished its growth during 
the oligarchic growth phase, so it is a stranded embryo that never participated in the giant impact phase \Jacobson and\ 
\Morbidelli\\20\^ . 

In the outer Solar System, the embryos need to grow to ~10 M 0 giant planet cores, so that gas accretion oc- 
1996|l. The core accretion scenario of dust growth, runaway growth and oligarchic growth is too 


curs ( Pollack et al. 


slow (for a review of the difficulties in this scenario, see Levison et al. 2010|l. Recently, an alternative solution has 
been proposed called pebble accretion {Ormel and Klahr 2010 Johansen and Lacerda\ |20T0{ Murray-Clay et al. 


\20\\\'^mbrechts and Johansen\\20\2\\Chambers 2014| l. While this pebble growth process was developed to solve a 

problem in the outer Solar System, it has strong implications for terrestrial planet formation. 

1.3 Pebble growth processes 

In the terrestrial protoplanetary disk, radial drift and the bouncing barrier frustrate growth at boulder sizes prevent¬ 
ing planetesimal formation. These barriers may be overcome by what may broadly referred to as gravo-turbulent 
growth (for a thorough review, stt \Johansen et a/^|2014| l, a possibly misleading nomenclature since turbulence isn’t 
always necessary. Gravo-turbulent mechanisms bring together such great particle concentrations that the collection of 
small particles collapses under self-gravity. By gathering enough dust or pebbles, a single body forms, which is signih- 
cantly larger than either of the barriers, skipping the boulder size regime altogether. The initial asteroid size-frequency 
distribution of the Main Belt may have been very shallow and dominated by asteroids larger than 100 km ^MorbidelU 


\et a/.|[2009] l, but this is not conclusive since growth from km sized asteroids can also match the current size-frequency 
distribution of the Main Belt {Weidenschilling\\20\ 1| |. However, there is also evidence from the Kuiper Belt that large 
equal mass binary systems with similar surface appearances were likely formed from a gravitationally collapsing cloud 
of particles that has too much angular momentum for a single body {Nesvorny et al. 2010|l. These clues suggest that 
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planetesimals appeared as 100-1000 km bodies. 

The importance of gas-particle interactions have a long history (e.g. Goldreich and War^ 1973t Weidenschilling 


\and CMzii||1993| l, but older ideas focused on creating a very thin, dynamically cold mid-plane layer of planetesimals. 
This disk could undergo a gravitational instability, but turbulent coupling with the gas and dynamical self-stirring 
makes this process difficult (Vifei£/ensc/!!7//ng anii Cmzz/ 1993 Youdin and Shu 2002|l. The more recently developed 


pebble mechanisms that allow growth beyond the boulder barriers fall into two categories. First, gas dynamics drive 
particles together via ‘passive’ concentration mechanisms. In regions of the disk with significant turbulence, eddies 
form at the smallest turbulent scales and these drive dust into the regions between the eddies ^Cuzzi et al.\ |2008[ 
\Chamber^ |2010| l. This mechanism creates smaller planetesimals, which may not be consistent with the evidence 
above. Larger vortices or radial velocity variations created by gap-clearing giant planets, phase transitions in the disk 
chemistry, or ionization levels of the disk, so called ‘dead zones’ create pressure bumps that collect particles (for 

2014|l. A pressure bump creates a zonal flow which slows 


a list of references see section 4.1.4 in Johansen et al. 


and possibly reverses the radial drift due to gas drag of pebbles ^Whipple\\\912\ . As mentioned, the pressure bumps 
are associated with specific structures in the disk, so if this mechanism is dominant, growth beyond the boulder 
barriers only occurs at unique locations. How numerous these locations are and where they are in the disk are open 
questions, but the answers so far appear to be few and typically they rest outside of the terrestrial planet formation 
region {Johansen et al. 2014| l. 

The second concentration mechanism is not associated with a pre-existing disk structure, instead the ‘streaming 
instability’ is a consequence of pebble and gas interactions. A single pebble causes almost no back-reaction on the gas 
dynamics, but a small clustering of pebbles causes the local gas to move with the pebbles at Keplerian orbital rates. 
Exterior pebbles catch up with this cluster via radial drift and increase the size of the back-reaction, further retarding 


the radial drift of the cluster and increasing the number of pebbles that catch up with the cluster {Johansen et al. 


20071. Eventually the cluster grows large enough to collapse due to self-gravity creating 100 to 1000 km directly from 


pebbles (for a thorough review of the material in the last two paragraphs, see \Johansen et al 2014| l. 

Pebble-pile planetesimals are an open topic of research ^Hopkins] |2014[ \Jansson and Johansen] |2014| l, but if this 
mechanism is important there are exciting repercussions. Eor instance, a fundamental unit of the Solar System is 
chondrule. Chondrules are small melt spherules found in the most abundant meteorite classes. If chondrules are the 
pebbles in these processes that would wonderfully explain their ubiquity, but currently most pebble concentration 
models require particles to be about an order or two of magnitude larger than a typical chondrule given standard disk 
models Johansen et a/. [[20071 1. Another repercussion of pebble concentration models are that the variations between 
planetesimal compositions may be a consequence of where in the disk it is most likely for the streaming instability to 
occur at a particular time. Due to radial drift, pebbles are constantly migrating from the exterior of the disk inwards. If 
there is an initial compositional gradient or a temporal gradient associated with dust sedimentation, then this gradient 
may be frozen in as a function of streaming instability location. Alternatively, pebbles may be the fragments from 
collisions between larger bodies { Kobayashi and Tanaka\^nW\\Chambers\\2Ql^ or fluffy aggregates created from icy 
dust grains ^Wada et al.\\2QQS]\Okuzumi et al 2012| l. It is not clear what the dominant pebble creation process is or, 
even, their appearances. 

Pebbles may also be important for the rapid growth of planetary embryos {Ormel and Klahr 2010[ Johansen 


\and Lacerda\ |2010[ \Murray-Clay et aL\ |201 1[ \Lambrechts and Johansen] |2012[ Chamber^ 2014| l. Pebbles are not 

entrained along streamlines in the nebular gas but neither are their motions independent of the gas. This has conse¬ 
quences for embryo accretion. As described above, planetesimals are accreted by embryos according to their mutual 
gravitational cross-sections, however this does not account for gas-solid interactions in the presence of a growing em¬ 
bryo. Pebbles can be accreted from a much larger cross-section than the gravitational focusing cross-section because 
as they react to the gravitational acceleration from the embryo, they experience gas drag, which reduces their velocity 
relative to the embryo and this increases the time at which the embryo’s gravity acts on them, and so on. 

In the outer Solar System, this process may allow the formation of the giant cores early enough to accrete gas 
envelopes \Lambrechts and Johanse n 2014 Lambrechts et al.\\20\A\ , although the cores could set-up a different per¬ 
sistent oligarchic growth regime ( Kretke andLevison 2014]l. In the inner Solar System, it also could play a significant 
role ( |Lev/5o« et al.\ 20141. If pebbles are a significant portion of the accreted mass, then embryos could become lay¬ 
ered with pebbles. If these pebbles are chondrules or clusters of chondrules, then these layers may be an explanation 
for the remnant magnetic fields and other evidence that certain chondrites were at the surface of larger bodies (this 
also assumes that these growing embryos have internally generated magnetic fields; |lTeA.!f et al.\ |2010| \Carporz^ 

201 l|l. Then again, these larger bodies could have formed directly from pebbles via 


et al. 


2011 Elkins-Tanton et al. 


a passive concentration or streaming instability mechanism. Most radically, if the streaming instability is inefficient. 
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then, perhaps, the five largest planetesimals grew via pebble accretion into the four terrestrial planets and Theia, the 
Moon-forming impactor—these are the only planet-sized bodies in the inner Solar System for which we have direct 
evidence of their existence. 

The consequences for pebble processes are revolutionizing planetesimal formation and the giant planets, but their 
repercussions for the terrestrial planets have not been fully explored. 

1.4 New pebble model from dust to embryo 

Assuming that pebble processes are the key to bypassing the boulder barriers, modified planet formation timeline 
stages emerge. Pre-existing grains or condensates settle into a dust structure centered on the mid-plane, and grow to 
pebble and boulder sizes from binary collisions. However, some as-yet-unknown fraction of pebbles are concentrated 
by the streaming instability and collapse directly into large 100—1000 km planetesimals. If this process is successful, 
then runaway growth occurs from this sea of planetesimals and planetary embryos emerge as oligarchs. If the popu¬ 
lation of planetesimals remains sparse, i.e. the streaming instability is rare, then embryo formation relies upon pebble 
accretion. Since the efficiency of pebble accretion is a sharp function of the size of the protoplanet, then planetesimals 
can be stranded at small sizes never experiencing rapid pebble accretion before the gas in the disk dissipates \Lam-\ 
\brechts and Johansen\\20\’^ . Interestingly, a bi-modal mass distribution between planetesimals and embryos seems 
inevitable regardless of the growth paradigm. In either paradigm, though, it’s unclear what the distribution of embryo 
masses are and how much mass is in the embryo population compared to the planetesimal population. 

This new model regarding the early stages of planet formation has strong implications for how possible initial 
compositional gradients are emplaced in the disk, for the timing of different accretion events relative to the decay 
of radioactive elements, etc. Thus, it requires much future study, but rather than waiting until these early stages 
are completely understood, studies of terrestrial planet formation have plunged ahead experimenting with different 
assumed initial distributions of planetesimals and embryos. 


2 Models of the giant impact phase of terrestrial planet formation 

Unlike earlier stages, which are often studied in local patches of the disk or in small annuli, the giant impact stage 
of terrestrial planet formation is usually modeled as a nearly global simulation. It is ‘nearly global’ because most 
N-body simulations simplify the treatment of the giant planets. They are often assumed to be fully formed, because 
they needed to accrete their gaseous envelopes during the lifetime of the gas disk, and often Uranus and Neptune are 
neglected since their influence on the terrestrial disk is weak (e.g. this was found to be so for the Grand Tack; \Walsh\ 

FSri[20TT] l. 

The terrestrial protoplanetary disk at this time is assumed to have a bi-modal mass distribution containing planetary 
embryos and planetesimals. Such a distribution is produced by both the ordered, runaway, and oligarchic growth 
scenario and a pebble-dominated scenario, and it is not clear which scenario is correct. The key parameter is the ratio 
of total mass in the embryo population to that of the total mass in the planetesimal population. It determines the amount 
of dynamical friction in the disk. Assuming ordered planetesimal growth, N-body simulations of the planetary embryo 
growth stage predicted a ratio near unity ^Kokubo and Ida\ 1998| l, but planetesimal grinding may significantly enhance 


this ratio over time ([LevAow et al. 20121 or perhaps, the planetesimals were only ever very small and collisional 


grinding and radial drift leaves behind a population of only Mars-sized embryos jKobayashi and Dauphas 20131. For 
pebble accretion, this ratio is unknown. Modem numerical N-body simulations have experimented with a large number 
of ratio choices {Jacobson and Morbidelli\\20\A\ or self-consistently model all the stages of planet formation ^Walsh\ 


and Levison\ 20141, and they are concluding that a high total embryo mass relative to the total planetesimal mass best 


reproduces the late accretion record on Earth (Jacobson et al. 2014| l. 

It is the orbits of Jupiter and Saturn that primarily distinguish the different models of the giant impact phase of 
terrestrial planet formation. The dynamical excitation of the giant planets is directly reflected in the dynamics of 
the terrestrial disk via secular resonances (particularly, and vq) and scattered embryos Raymond et a/.|[2009| l. We 
discuss three classical models, in each the terrestrial disk extends from an inner edge to an outer edge carved by Jupiter. 
In the classical scenarios, the giant planets do not migrate, but they have significantly different orbits in each model 
due to different assumptions about past and future evolution of the Solar System. Each assumption has consequences 
for matching constraints such as the late heavy bombardment. Truncated disk models assume that either the terrestrial 
disk does not extend out to the giant planets or that the migration of the giant planets truncates the disk. This Grand 
Tack model is the most successful and is described in great detail. 
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Figure 1; The mass and semi-major axis distribution of planets formed from the eccentric (empty squares; O’Brien 


et al. 

2006 

et al. 

2009 




aymond et al. 2009 Fischer and Ciesl^\2Q14\, circular (empty circles; [O’fiWew et aL\\2QQ6\\Raymond 


20141, and extra-eccentric (empty triangles; Raymond et al. |2009|l Jupiter and Saturn 


^ischer and Ciesla 

models and the Grand Tack model (solid points; Walsh et al.\ 201 1[ O’Brien et al. 2014[ \facobson et al. 2014 


Jacobson and Morbidelli 2014|l. Planets within the upper gray rectangle are considered Earth-like and those within 


the lower gray rectangle are considered Mars-like. 


Comparing terrestrial planet systems Numerical simulations never exactly reproduce the terrestrial planets and 
cross-comparing the simulated planets with the true distribution of masses and orbits one-by-one is a complicated 
task. While, we can directly compare the orbits and masses of the planets in each system as shown in Figure 1 
evaluating which model best reproduces the Solar System is a matter of statistics and gross metrics. To help with this 
process, Chambers (2001|l started using the following sophisticated quantities. 


The angular momentum dehcit Sd of the inner terrestrial planets measures their dynamical excitation [Laskar 


|1997[ ). It is dehned as the difference between the sums of the angular momentum of the corresponding circular and 
in-plane orbits and the actual orbit normalized by the sum of the angular momentum of the corresponding circular and 


in-plane orbits [Chambers 20011: 


Sd = 






— \/l — Cj cos ij 




( 1 ) 


where ruj is the mass, aj is the semi-major axis, Cj is the eccentricity and ij is the inclination of the jth planet. The 
inner planets do strongly perturb one another exchanging angular momentum, but exchanges with the outer planets are 


limited so the quantity Sd is conserved within a factor of two [Laskar 19971. The current value for the Solar System 
is Sd = 0.0018 and is marked by a dashed horizontal line in Figure 2. If the Solar System undergoes a giant planet 
instability as supposed in the Nice model, then the angular momentum dehcit after terrestrial planet formation needs 
to be between 10% (0.00018) and 70% (0.00126) of the dehcit after the instability (shown in Figure 2 as the gray 
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Figure 2; The angular momentum deficit Sd and concentration statistic Sc for terrestrial planet systems formed from 
the eccentric (empty squares; O’Brien et al.||2006 Raymond et al. 2009 Fischer and Ciesla 2014) l, circular (empty 
circles; lO’Bnen et al. 2006[ Raymond et al.\ 2009 \Fischer and Ciesla 2014|l, and extra-eccentric (empty triangles; 


Raymo^ et a/.||2009|l Jupiter and Saturn models and the Grand Tack model (solid points et ar\\20\\\\b’Brien\ 


et al. 2014 Jacobson et al.\^0\A\\Jacobson and MorbiddU\\20\^. The horizontal dashed line is the angular mo¬ 


mentum deficit of the current terrestrial planet system Sd = 0.0018 and so the target value for the eccentric and 
extra-eccentric Jupiter and Saturn models. The gray rectangle marks the region with 10-70% of the current angular 
momentum deficit ^rasser et al.\ 2013[ ), and so is the target for models like the Grand Tack and circular Jupiter and 
Saturn that include a giant planet instability. The vertical solid line indicates the current terrestrial planet system’s 
concentration statistic Sc = 89.9. 


rectangle; Brasser et al. 2013| l. 

The other useful but unusual quantity is the concentration statistic Sc. The terrestrial planetary system of the Solar 
System is interesting because its mass is almost entirely in Earth and Venus with only a little bit in Mercury and Mars. 
The concentration statistic attempts, albeit degenerately, to capture this mass-orbit distribution with a single number. 
It is defined ( | Chamber^ |200 1 [ ) : 


Sr = max 


Jljruj (logio («/aj))' 


( 2 ) 


where like before roj is the mass and Oj is the semi-major axis of the jth planet. Sc is the maximum of the bracketed 
quantity as a is varied. The terrestrial planets in the Solar System have Sc = 89.9, and this value is marked by a solid 
vertical line in Figure 2. 


2.1 Classical models 


In the classical models, the terrestrial disk of planetesimals and embryos has an outer boundary set by Jupiter, since 
most orbits exterior to the inner 3:2 mean-motion resonance with Jupiter are unstable. However, the inner boundary 
is more mysterious, and there is no conclusive theory regarding why the Solar System doesn’t have planets interior to 
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Mercury. Existing theories include high collision velocities that fragment and grind away the planetesimal population 
before embryos can form ( Chamber^ 2001| l, gas driven migration of embryos into the Sun ^da and Lin\ |200^ , an 
inner edge to the gas disk, or a silicate evaporation front that doesn’t allow dust or pebbles to exist interior to form 
planetesimals and then embryos. Most N-body simulations place the inner edge at 0.3-0.7 AU to avoid forming planets 
interior to Mercury. Forming a Mercury-sized planet in the correct orbit with a satisfactory explanation for its high 
metal content is an open problem. 

Eccentric Jupiter and Saturn, current orbits The first giant impact phase models considered disks that contained 
about half the mass in the planetesimal population and half in the embryo population, which is consistent with the 
mass ratio found when oligarchic disks go unstable ^okubo and Ida\ |1998[ ). Due to dynamical friction, these bodies 
start on low eccentricity and low inclination orbits. Importantly, Jupiter and Saturn are assumed to have formed on 
their current orbits, and the N-body simulations begin after the gas has been removed from the disk. 

This model is remarkably successful, it provided confidence that these models in general are on the right path. 
Early numerical simulations that only included a limited number of embryos naturally formed about four planets near 
the correct semi-major axes of the terrestrial planets (see Figure 1), giant impacts similar to the Moon-forming impact, 
and formation timescale of tens of millions of years for Earth-like planets y^^gnor et al 1999 Chamber^ 2001) 1. 
When enough planetesimals are included in the simulations, dynamical friction removes excess energy and angular 
momentum from the orbits of the growing planets transferring it to the planetesimal population and the simulations 
reproduce the low eccentricities and low inclinations of the terrestrial planets ^O’Brien et aL\ |2006[ \Monshima et al.\ 
|2008 |). Some but not all terrestrial planet systems created by the eccentric Jupiter and Saturn model have angular 
momentum deficits consistent with the Solar System \0 ’Brien et al. 200b\^aymond etW. 2009) Fischer and Ciesla 


20141; since the eccentric Jupiter and Saturn model is inconsistent with a Nice model-like giant planet instability, the 


most consistent angular momentum deficits are those closest to that of the current terrestrial planets (the dashed line 
in Figure 2). 

This model is haphazardly successful at delivering water to Earth-like planets from the outer asteroid belt via the 
vq secular resonance y\4orbidelli et al.\ |2000| {Raymond et al.\ |2004[ {O’Brien et al.\ |2006[ {Raymond et al.{ |2009| ). If 
most of Earth’s water is delivered from the planetesimal carbonaceous chondrite parent bodies, then the D/H ratio of 
Earth’s water is naturally explained (see Morbidelli et al.\ 2000 for greater detail and a discussion of the less likely 


alternatives). Inspired by the data from primitive meteorites, a simple model for the water content of planetesimals as 
a function of radius is that the water mass fraction is 10“^% interior of 2 AU, 10“^% between 2 and 2.5 AU, and 5% 
exterior to 2.5 AU ^Raymond et al.[ [2004) 1. This model is used for all of the classical scenario simulations in Figure 
3, but not for the Grand Tack simulations since that model assumes that the carbonaceous chondrite parent bodies 
originate from exterior to Jupiter and Saturn. 

There are four outstanding issues with the eccentric Jupiter and Saturn model. First, this model systematically 
creates a planet at the position of Mars with a mass between 2 and 12 times that of Mars. As shown in Figure 4, there 
are a few exceptions to this rule and this directly leads to low values of the concentration statistic as shown in Figure 2. 
These were discovered by simulating a much larger number of simulations than had historically been studied before, 
emphasizing the need for this approach y^ischer and Ciesla 20\A\ . Of the 62 numerical simulations, 6 (^10%) have 
Mars-like planets exterior of all Earth-like planets in the system. If Mars is an outlier of the eccentric Jupiter and 
Saturn model, then this is a challenge to the Copernican principle—that the Solar System is a usual rather than special 
outcome of planet formation processes. Before such a conclusion should be drawn, a more thorough investigation into 
what happened differently in these 6 simulations needs to be completed. 

Raymond et al. ( 2009)1 demonstrated that imperfect accretion via high velocity impacts cannot be responsible for 


the difference between the masses of the planets formed at Mars location and the mass of Mars. This was confirmed 
by Chambers (20131, who found that while imperfect accretion makes smaller planets overall including at Mars’ 
location, it tends to decrease the mass ratio between Earth and Mars rather than increasing it relative to the perfect 
accretion simulations. Also, unlike Mercury, Mars appears inconsistent with this type of impact-driven mass loss. 

Second, it does not explain how Jupiter and Saturn obtained their current eccentricities. Interactions with the nebu¬ 
lar gas disk should leave the gas giants on quasi-circular orbits ( Kley 2000) Bryden etal. \ 2000) Snellgrove et al. 2001 


Papaloizou 2003) Crida et ^)2008)l. Also, giant planet migration is likely to move Jupiter and Saturn into a mean 


motion resonant configuration, and not into the near 5:2 mean motion configuration they currently occupy y^apaloiz^ 
{et al.{^b0\\\Kley and Dirksen{{200^\Bitsch et al.))2013) ). 

Third, this model is inconsistent with the Nice model, and so does not directly account for the late heavy bombard¬ 
ment. As discussed in Section o alternative explanations for the late heavy bombardment exist although none are 
as complete or successful as the Nice model. Only if an alternative is accepted, can the eccentric Jupiter and Saturn 
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Figure 3: The water mass fraction for each Earth-like planet formed from the eccentric (empty squares; Raymond 
\et al.\ |2009[ ^ischer and Ciesl^ |2014 1, circular (empty circles; [RayOTonJ et al.\ | 2009[ ^ischer and Ciesl^ |2014 1, and 
extra-eccentric (empty triangles^J^flymow^^^^ 2009 i Jupiter and Saturn models and the Grand Tack model (solid 


points; Walsh et g/T] 201 1[ O’Brien et al. 2014 Jacobson et al.\ 2014 Jacobson andMorbidelli 2014[ ) as a function 
of the mass of the planet. The water mass fraction is the water mass delivered to the planet assuming perfect accretion 
and either the water model of \Raymond et aZ.| ( |2004| l for the non-Grand Tack simulations or the model Qf \0’Brien et al.\ 
( 2014| l for the Grand Tack simulations. Earth-like planets from the numerical simulations have masses between 0.5 
and 2 M 0 and orbits between the current orbits of Mercury and Mars as shown in Eigure[^ The upper dashed line is 
a liberal estimate of Earth’s water content ( [Marty 2012 1 , the gray region is a more probable estimate \Lecuyer et al.\ 
|1998| ), and the lower dashed line is a minimum { Lecuyer et aZ.|[T998] l, i.e. an estimate of the known water reservoirs. 
The water content of Earth’s mantle reservoir is unknown. 


model be consistent with the evidence for the late heavy bombardment. 

Eourth, the eccentric Jupiter and Saturn model is internally inconsistent since interactions with planetesimals and 
embryos in the terrestrial disk tend to circularize the giant planets ( Chambers 2001) l. 

Extra-eccentric Jupiter and Saturn The extra-eccentric Jupiter and Saturn model was devised by Chambers and 
|Ca.!f.!few| ( |2002j l to alleviate the fourth outstanding issue of the previous model. The giant planet are assumed to have 
formed at their current semi-major axes but are on more excited orbits. This dynamical excitation is transferred through 
both scattering events and secular resonances, principally the 1/5 and the to the terrestrial disk, and this creates an¬ 


gular momentum deficits much higher than those possessed by current terrestrial planets in the Solar System {Raymond 
et al. \ 2009 j l, as shown in Figure 2. 


The rapid depletion of the outer asteroid belt typically leads to a very dry Earth inconsistent with current water 
fraction estimates as shown in Eigure 3. If the gas disk disperses such that Mars mass embryos in the outer asteroid belt 
region, assumed to be water-rich, are swept up in evolving secular resonances with Jupiter, then they may be incorpo¬ 
rated into Earth delivering the needed water ( Thommes et al. 2008| Raymond et al.\ 20091. Alternatively, water may 


be delivered by mechanisms listed by Raymond et al. (20091 such as water adsorption of small silicate grains {Mu- 


\ralidharan et al.\ |2008jl, comet impacts ( Owen and Bar-Nun\\\995\, or oxidation of a H-rich atmosphere \Ikoma and 
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Figure 4: Lines connect the final planets of the good terrestrial system analogs. Otherwise the individual graphs are 
very similar to FigureFor the classical models, good terrestrial system analogs have at least one Mars-like planet. 
For the Grand Tack model, good terrestrial system analogs have at least one Mars-like planet and both a Venus-like 
planet and Earth-like planet with no planets in between. 


\Genda\\2(X)(i) . 

Sweeping secular resonances combined with tidal gas drag can effectively deplete the asteroid belt by driving the 
inward migration of planetary embryos ^agasawa et al.\ | 2005[ ) and reproducing the small mass of Mars ( Thomme^ 
|2008|l, but this requires a rather specifically tuned gas disk, particularly the gas removal rate ^Thommes et al. 


et al. 


20081. Sweeping resonances due to changing precession rates as the gas is dispersed is not necessary, fixed frequency 


resonant and secular perturbations do clear out the asteroid belt and Mars-region occasionally creating a small Mars, 
as shown in Figure 4, however there are often planetary embryos remaining in the asteroid belt Raymond et al 2009) . 
Thus, the extra-eccentric Jupiter and Saturn scenario weakly alleviates the first outstanding issue identified above, but 
many final planets in the Mars region are still too large and since the excitation retards the growth of Earth-like planets, 
they tend to be too numerous and too small. This leads to unchanged concentration statistics as shown in Eigure 2. 

The extra-eccentric Jupiter and Saturn scenario accentuates the second problem identified with the current orbit 
scenario; the assumption of high initial eccentricities for the giant planets is inconsistent with models of their growth 
in a gas disk ^Morbidelli and Crida\\2007\\Pierens and Raymond\^011\\Bitsch et aZ.|[2013| l. Similar to the eccentric 
Jupiter and Saturn scenario, this model is not coherent with the Nice model. 

Einally j^Qymowt/ et al.\ ( 2009| l identifies a number of other explanations for extra-eccentric giant planets such as 
an early giant planet instability leaving the giant planets on very eccentric orbits and the corresponding problems with 
these models. 
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Circular Jupiter and Saturn, pre-Nice 2.0 model Numerical simulations consistent with the Nice model automat¬ 
ically solve all of the latter three outstanding issues of the eccentric Jupiter and Saturn model {Raymond et al. 2009| l. 
Giant planets growing and migrating in a gas disk circularize and enter resonances, and it is natural for Saturn and 
Jupiter to enter into a 3:2 mean motion resonance near 5.4 AU and 7.3 AU ^Morbidelli et a/.|[2007[ ). This configuration 
of the giant planets is primed for a giant planet instability at ^400 My due to interactions with a massive outer disk of 
which the Kuiper Belt is a remnant {Levison et al. 201 1| |. 

This model has many of the same successes as the eccentric scenario: correct number of planets including Earth¬ 


like planets near 1 AU as shown in Figure 1, giant impacts similar to the Moon-forming impact { Raymond et al.\ 
|2009| ), and growth timescales of tens of millions of years Raymond et al.\ |2009[ ^ischer and Ciesla J20M|. Since 
the giant planets are on circular orbits, excitation due to secular resonances is reduced ( Raymond et al. |2009 ^ischer 
and CiesJ^ 2014|l but since the giant planet instability in the Nice model will dynamically excite those planets, the 


target angular momentum dehcit value is 10-70% lower than for the eccentric Jupiter and Saturn model {Brasser et al. 


|2013| ). Such systems are successfully created as shown in Figure 2. Another success is the delivery of water to Earth 
from the asteroid belt. Using the same model as before {Raymond et al. 2004| l, enough water from the outer asteroid 
belt is always delivered to Earth as shown in Figure 3 Raymond et al.\ 2009 Fischer and Ciesla\\20\A) . 

Notably though, the small Mars problem is accentuated, as shown in Figure 4. Planets at the location of Mars 
are almost all between 2 and 20 times too massive. Only 1 (~2%) of 62 systems has a Mars-like planet without also 
containing a third near-Earth mass planet. This reveals itself also as very low concentration statistics, as shown in 
Figure 2. 

While this model is consistent with the Nice model, it is only partially consistent with evolution in a gas disk. The 
circularization and migration of the giant planets into mean motion resonances in the circular Jupiter and Saturn model 
is assumed to have no effect on the evolution of the terrestrial disk. The Grand Tack model is the natural consequence 
of breaking this assumption. 

2.2 Truncated disk models 

As mentioned, one of the outstanding problems with the classical models of terrestrial planet formation is the size of 
the planet formed near 1.5 AU. As shown in Figure 1, these models systematically produce planets at the location 
of Mars that are about 5-10 times too massive. Raymond et a/.| ( 2009| l considered a broad range of models including 
those discussed above and a few others, but none consistently solved this problem and most introduced other issues 


as well. This problem is an old problem that has been recognized at least as far back as Wetherill (1991 1. Walsh and 
\MorbidelJl\^20l l| l considered planet migration and the resultant sweeping secular resonances as a way to deplete the 
region around 1.5 AU of solid material early enough to frustrate the accretion of Mars-analogs. This failed as the 
migration of the giant planets would happen too late compared to the accretion of Mars-like planets, leaving them at 
similar mass as in classical scenarios. 

Hansen (2009|l proposed a dramatic solution to this problem demanding wholesale changes to the surface density 


of the solid material in the disk. They found that if the initial disk is truncated at 1 AU then the mass distribution of 
the terrestrial planets could be recovered including the Earth/Mars mass ratio. The systematic success to the Hansen 
model stems not just from the disk truncation but the unstable dynamical over-packing of the 400 equal mass embryos, 
which begin the simulation with overlapping mutual Hill radii ^O’Brien et q/.[| 2OT4) |. These initial conditions lead to 
the rapid scattering of some embryos out of the narrow annulus. The bodies scattered out beyond 1 AU were typically 
stabilized due to interactions with other scattered bodies and isolated beyond the edge of the disk. This allowed these 
bodies to stop their accretion early and remain at a low mass, resulting in a short accretion time relative to the planets 
that accrete within the annulus. This is similar physics to that observed in previous works that had disk edges around 
1 or 1.5 AU due to computational requirements ^gnor et al.\ 1999\ Chambers 200 1[). 


Jin et al. (20081 proposed that severe surface density depletions in the proto-planetary disk cause a truncation. 


Since different angular momentum transport mechanisms are active at different radii in the disk, [fin et al. (20081 
found a non-monotonic gas surface density for the inner regions of the disk and a surface density minimum at about 
1.5 A\J. \Izidoro et al. \\20\A \ performed N-body experiments and deduced that to best fit Mars, a ^ 75% surface 
depletion between 1.3 and 2 AU is needed, however the prediction by Jin et al. (20081 appear to be more consistent 
with ~ 25% and over a narrower annulus. Unlike the hypothesis proposed in Jin et al. ( 2008| l, Izidoro et al. ( 2014| l hnd 
that Mars is a stranded embryo scattered out from the interior region of higher surface density, similar to the Grand 
Tack. Dissimilar to the Grand Tack and like the classical scenario, the asteroid belt is grown from material already in 
the low surface density region in the|y/» et a/.|(|2008^ and\Izidoro et a/.|(|2014|l models. 
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‘Grand Tack’, migrating Jupiter and Saturn A narrow over-packed disk with a truncation of solid material at 
1 AU produces ideal outcomes for the Earth/Mars mass ratio { Hansen\ 20091, but why such a disk existed needs an 
explanation. However, a powerful means to re-distribute solid material in the early Solar System are the giant planets. 
The giant planets must have had a relatively short accretion timescale compared to that of the terrestrial planets (see 
Section [13 and so existed to move around terrestrial building block material. Hydrodynamic models and data from 
extra-solar planetary systems suggests that planets the size of Jupiter can carve an annular gap in a gaseous nebula 
resulting in their inward migration due to the viscous evolution of the disk ( Ward\ 1980[ Lin and Papaloizou 1986) 1. 
This migration, type-II migration, can lead to the inward migration of a Jupiter-mass planet on 100,000 year timescales. 

A series of works found that the presence of a second massive planet, a Saturn mass planet, can halt and even 
reverse the inward migration of a Jupiter mass planet ^Masset and Snellgrove 200]| Pierens and Nelson 2008 Pierens 


\et al.\ |2014[ \Crida and Morbidelli] |2007[ \Pierens and Raymond |201 1[ D’Angelo and Marzari\ 2012| l. Saturn is found 
to migrate very close to Jupiter, reaching an exterior 2:3 mean motion resonance with Jupiter, and this proximity causes 
their gaps to overlap. The net effect is that Saturn’s presence allows gas to pass through Jupiter’s gap, stopping the 
inward type II migration. The lower surface density of gas behind Jupiter, owing to the presence of Saturn’s gap, then 
allows for the inner disk to provide a much stronger torque on Jupiter’s orbit pushing both planets outward. 

The net effect of this inward-then-outward migration scenario is that it provides a possible mechanism to create 


a truncated disk of solids. Walsh et al. (2011 1 found that if Jupiter reversed its inward migration at 1.5 AU, i.e. 
‘tacked’, then the disk of solids would be truncated at 1 AU, however the violent mechanism causing the truncation 
and over-packing the disk leaves a more excited initial state than that previously explored by )HQ«5e7i| ( |2009 |l. Similarly, 
while ^Iansen\ ([2009]l started with 400 similar-size bodies in the annulus, l\\t\Walsh et a/.|(|201 11 model, dubbed ‘the 


Grand Tack’, started with a bi-modal size distribution of embryos and planetesimals. Despite these differing initial 
conditions, the terrestrial planets produced were similar including the high Earth/Mars mass ratio as shown in Eigure 
1 . 

The Grand Tack creates an asteroid belt very differently than the classical models since Jupiter migrates through 
the asteroid belt-twice. Rather than depleting a pre-existing belt via secular resonances and scattering events, the 
asteroid belt in the Grand Tack is created by the inward scattering of material by the giant planets as they migrate 
outwards, but not all this material was originally from the outer Solar System. Some of it had been scattered into 


the outer Solar System from the inner Solar System originally. Walsh et al. (2011 1 hypothesize that material from the 
outer Solar System corresponds to the broad C-complex of asteroids associated with the carbonaceous chondrites and 
the material that originated in the inner Solar System but was scattered out and then back in corresponds to the broad 
S-complex of asteroids associated with the ordinary chondrites. The classical models explain the observed, albeit 
messy, gradient from S-complex to C-complex in the asteroid belt [DeMeo and Carry 2013) 2014| l as a local property 
of the disk. However, the Grand Tack asserts that this transition occurs over a much larger radial section of the disk 
(interior and exterior to the giant planets), but the migration of the giant planets naturally re-creates this gradient in the 
asteroid belt ^Walsh et al.||2011]|2012| ). 

This hypothesis has strong implications for delivery of water to Earth since the D/H ratios of carbonaceous chon¬ 
drites best match Earth’s water ^Morbidelli et al. \ 20001. These outer planetesimals account for a total of 1-3% of the 
total accreted mass for each of the terrestrial planets ( Walsh et al. 2011) \O^Brien ef oT 2014| l. As shown in Eigure 3, 
they deliver enough water to Earth assuming that they contain ^10% water by mass. Moreover, unlike the classical 
models, the mass of C-complex material delivered to the planets is independently calibrated because the initial ratio 
of C-complex to S-complex material in the protoplanetary disk is known from the ratio that survives in the asteroid 
belt {DeMeo and Carry |2013 2014| l. 

The Grand Tack solves each of the outstanding problems identified with the eccentric Jupiter and Saturn model. 
The giant planets and the terrestrial disk evolve self-consistently in a gas disk and the Nice model is naturally incor¬ 
porated. Venus-like, Earth-like and Mars-like planets are all natural outcomes and many systems have an architecture 
similar to the Solar System as shown in Eigure 4. This is why the Grand Tack alone matches both the concentration 
statistic and the necessary angular momentum deficit, pre-Nice model in this case, as shown in Eigure 2. 

There are still some open issues regarding the Grand Tack such as the formation of Mercury, which may just be a 
result of numerical resolution, and some discussed in \Raymond and Morbidelli\^2014\ . 

3 Earth in the Grand Tack model 

The Grand Tack model is best able to reproduce the mass-orbit distribution of the terrestrial planets and is the most 


consistent with the history of the Solar System. This includes reproducing the rapid growth of Mars ( Dauphas and 


Pourmand 2011) Nimmo and Klein^ 2007)l by acknowledging that Mars is a stranded embryo {Jacobson and Mor- 
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Figure 5; The relative late accreted mass for each Earth-like planet formed from the eccentric (empty squares; O’Brien 


et al. 


et al. 


2006[ \Raymond et al. 2009 


Fischer and Ciesla 


2009 Fischer and Ciesla 2014| l, circular (emp ty circles; [O’fi Wew et al. 2006[ Raymond 
20091 [fFfc/ier Ciesla 20141, and extra-eccentric (empty triangles; Raymond et al. |2009| l Jupiter and Saturn 

models and the Grand Tack model (solid points; \Walsh et al. \ |201 1[ \0’Brien et al. 2014[ [ygcofeow et al.\ |2014[ 
\Jacobson and Morbidelli\ |2014| l as a function of the time of the last giant impact. The relative late accreted mass is 
the mass accreted after the last giant (embryo) impact divided by the final mass of the planet. Earth-like planets from 
the numerical simulations have masses between 0.5 and 2 M 0 and orbits between the current orbits of Mercury and 
Mars as shown in Eigure 


bidelli\ 2014| l. If Mars is an embryo then we have some information regarding the size of planetary embryos, and 
if planetesimals are born big, ^100-1000 km, as the new story regarding pebble formation suggests, then the only 
remaining significant free parameter describing the bi-modal mass distribution is the ratio of the total mass of the em¬ 
bryo population to the total mass of the planetesimal population. This ratio has direct consequences for the dynamical 
friction present in the terrestrial disk: more planetesimal mass means more dynamical friction [Ida andMakino 1993|l. 


Jacobson and Morbidelli (20141 showed that increased dynamical friction also leads to an earlier last giant impact 


since lower relative velocities accelerate growth. Thus, the later the last giant impact, which is also the Moon-forming 
impact, then the embryo population must be more massive relative to the planetesimal population. \ Jacobson et al.\ 
( |2014| l showed there is a strong correlation between the time of the last giant impact and the late accreted mass, 
which is the mass accreted by Earth-like planets after the last giant impact. This makes sense because the mass of 
planetesimals in the disk decays with time, so the amount of planetesimal mass accreted after a later time is smaller 
than after an earlier time. This correlation for both the Grand Tack and the classical models is shown in Eigure 5. 
Classical models always date the Moon-forming impact later than the Grand Tack. 

The geologic record in Earth’s mantle provides an estimate of the late accreted mass. From the highly siderophile 
element record on Earth (Becker ef aZ. 2006 Walker 20091, it appears that ~5x 10“^ Earth masses, Mq, of material 
was delivered to Earth after the Moon-forming event (many lines of evidence restrict this late accreted mass to below 
0.01 M 0 ; Jacobson et al. 2014 also see Morbidelli and Wood in this monograph). Using the correlation between the 
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late accreted mass and the time of the last giant impact, this late accreted mass estimate dates the Moon-forming impact 
to ^95 My {Jacobson et al. 2014| l. Such a late last giant impact requires that most of the mass in the terrestrial disk 
be in the embryo population rather than the planetesimal population when Jupiter’s migration interrupts the oligarchic 
growth phase \Jacobson and MorbidelU 2014jt^This is consistent with models of the oligarchic growth phase that 
include planetesimal grinding (Levison et al. 2 Q\2\ as well as a pebble model that rarely creates planetesimals via the 
streaming instability but efficiently turns planetesimals into embryos via pebble accretion. 

3.1 The growth of Earth 

We now delve deep into the details of the growth of Earth in the Grand Tack scenario. In the prior sections, we 
were exploring the overall differences between the terrestrial planet formation models. Here, we focus on specific 
simulations that produce terrestrial planet systems most like the Solar System. In order to do that we make stricter 
requirements. We require that a Venus-like planet with mass between 0.4075 and 1.63 M 0 be directly interior to 
an Earth-like planet with a mass between 0.5 and 2 Mq with at least one Mars-like planet exterior to both of those 
planets with a mass between 0.0535 and 0.114 Mq. Eurthermore, Earth-like planet most have a late accreted mass 
between 0.001 and 0.01 Mq. These stringent requirements leave only 10 planetary systems out of the 203 Grand 
Tack numerical simulations from \Walsh et al.\ ( |2011[ l, \0’Brien et al.\ ( |2014| l, \Jacobson et al.\ ( |2014| l, and \lacobson\ 
and MorbidelU ( 2014| l and all are from the latter two publications since they include the only simulations with a high 
enough total embryo to planetesimal mass ratio. Three of the ten have two Mars-like planets, but we include them for 
better statistics. 

The growth of these Earth-like planets are shown in Eigure 6. They all follow a very similar trajectory and grow 
quickly at first reaching 50% of their growth in Tso^lO My, however their growth is strongly non-linear. In the past, 
many have approximated this non-linear growth with an exponential function oc exp (f/r). The best hts to these 
simulations have exponential timescales t~3-10 My, but we find exponential fits to be inappropriate for the Grand 
Tack. Eigure 6 shows exponential hts to the growth prohle of each Earth-like planet, and then shows the difference 
between the actual growth trajectory and those hts. There is a clear underestimate of the growth rate at early times. A 
much better model is the Weibull model c>c exp [t/a)^. The best hts to these simulations are shown in Eigure 6 with 
a~10 My and /3~0.5. This model has no systematic over- or under-estimate. However, even this model can be off by 
20 %. 

3.2 Composition of Earth and the other terrestrial planets 

In the Grand Tack, the inner Solar System is thoroughly radially mixed by the inward migration of Jupiter (|0 ’Brien 


et al. 20141. This naturally produces Venus-like and Earth-like planets that are compositionally very similar. Both 


planets contain multiple embryos and many planetesimals from throughout the disk. Eigure 7 shows their compositions 
as a cumulative function of semi-major axis. There are three different compositional regions made clear. Eirst, ^75% 
of the mass originates from 0.7 AU, the inner edge of the terrestrial disk, to ^1.1 AU, the location of the 3:2 inner 
mean motion resonance with the location of Jupiter’s ‘tack’ at 1.5 AU. The mass from within this region is remarkably 
linearly representative. About ^22-25% of each planet’s mass originates from exterior of 1.1 AU to the outer edge of 
the terrestrial disk at about 4 AU. This mass is not delivered linearly with semi-major axis, but with slightly more mass 
delivered from the inner section than the outer. This is because material scattered by Jupiter further out in the disk was 
less likely to be scattered all the way into the 1 AU region where the proto-Earth and proto-Venus were growing. The 
hnal ~l-3% of mass comes from the outer disk exterior to the giant planets (O’Brien et al. 2014] l. 

The feeding zones of the proto-Earth before the last giant (Moon-forming) impact is shown for each simulation in 
Eigure 8. It very much resembles the hnal Earth-like planet, which makes sense because Theia, the last giant impactor, 
and the late accreted mass, which are also shown in Eigure 8, contribute only a small amount of mass. The late accreted 
mass appears to sample the entire inner disk. Unlike the hnal Earth-like planet, the composition of the late accreted 
mass doesn’t rehect the location of Jupiter’s tack. 

Theia is usually a stranded embryo like Mars, although sometimes Theia is much larger and looks very similar 
to the proto-Earth. The feeding zones of Mars-like planets are also shown in Eigure 8 and these reveal that while 
Mars-like planets often come from the outer disk, they do not always. Similarly, while Theia often comes from the 
inner terrestrial disk near 1 AU, it can come from the outer terrestrial disk. Some Mars-like planet compositions are 
composed of different embryos and so are unlikely to match the Hf-W growth constraints. 
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4 Conclusion and discussion 


The growth of Earth cannot be isolated from the formation of the rest of the Solar System. While the number of 
constraints regarding the history of the Solar System are few, they are powerful. The mass and orbit distribution of 
the terrestrial planets and the evidence for the late heavy bombardment rule out the classical terrestrial plant formation 
models. Trunfcated disk models reproduce the mass and orbit distribution of the terrestrial planets particularly the 
Earth/Mars mass ratio. Of these models, only the Grand Tack is completely self-consistent and consistent with a late 
heavy bombardment via the Nice model. 

The growth of Earth in the Grand Tack model follows a Weibull accretion curve to within 20%. During its growth, 
it samples ~75% of its mass linearly from the inner terrestrial disk edge to the location of the 3:2 inner mean motion 
resonance with the location of Jupiter’s tack. The remainder of its mass comes predominately from the rest of the 
terrestrial disk but biased towards the regions exterior but closes to 1.1 AU. Plenty of water is delivered to Earth via 
C-complex asteroids from the outer Solar System. The total mass of these asteroids delivered into the inner Solar 
System is calibrated by the ratio of S-complex to C-complex asteroids in the Main Belt. 

Erom the highly siderophile record on Earth, the age of the Moon is estimated to be ^95 My ^Jacobson et al. 


2014[). Since this age estimate is high, the ratio of the total mass in the embryo population to the total mass in the 


planetesimal population must also be high {Jacobson and Morbidelli\l20l4). Since Mars is a stranded embryo {Nimmo 


and Kleine 


Morbidelli 


2007 1 \Dauphas and Pourmand 201 1||, the mass of individual embryos is also known {Jacobson and 


2014|l. These two facts strongly constrain the structure of the terrestrial disk at the time of Jupiter’s inward 


migration. While this could be explained by planetesimal grinding, it could also be the result of pebble processes. 

Pebble processes are likely the key to growing past the boulder barrier and creating planetesimals via the streaming 
instability. They also likely contribute a significant amount of mass to embryos via pebble accretion. These processes 
are likely to leave unique compositional gradients and signatures in the disk. Combining these gradients with models 
of the composition of Earth and observations of Earth’s composition will be a powerful test of these hypotheses. 
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Figure 6; Growth curve analysis for Earth-like planets. Upper left graph shows the growth of Earth-like planets 
as described in Section HU and the left middle and bottom graphs are hts using exponential and Weibull models, 
respectively. The upper center plot shows the time to reach 50% of the hnal mass as a function of the exponential ht 
timescale and the upper right plot shows the a and /3 of each of the Weibull hts. The right middle and bottom graphs 
are the the difference between the growth curve and the exponential and Weibull hts, respectively. 
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Figure 7; The composition curves for Venus-like and Earth-like planets. These curves correspond to the terrestrial 
planetary systems described in Section 3.1 The composition curves show the cumulative mass in the final planet as a 
function of semi-major axis. Some mass is accreted from exterior of 3 AU. 
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Figure 8; The composition curves for the Theia-like, Mars-like, and proto-Earth-like planets and the late accreted mass 
composition. These curves correspond to the terrestrial planetary systems described in Section 3.1 The composition 
curves show the cumulative mass in the final planet as a function of semi-major axis. Some mass is accreted from 
exterior of 3 AU. 


18 
























































































































References 


Agnor, C. B., R. M. Canup, and H. F. Levison (1999), On the Character and Consequences of Large Impacts in the 
Late Stage of Terrestrial Planet Formation, ICARUS, 142(1), 219-237. 

Alexander, C. M. O. D., J. N. Grossman, D. S. Ebel, and F. J. Ciesla (2008), The Formation Conditions of Chondrules 
and Chondrites, Science, 320(5), 1617. 

Allegre, C. J., G. Manhes, and C. Gopel (2008), The major differentiation of the Earth at ^ij4.45 Ga, Earth and 
Planetary Science Letters, 267(1-2), 386-398. 

Armitage, P. J. (2014), Lecture notes on the formation and early evolution of planetary systems, pp. 1-71. 

Becker, H., M. E Horan, R. J. Walker, S. Gao, J. P. Lorand, and R. L. Rudnick (2006), Highly siderophile element 
composition of the Earth’s primitive upper mantle: Constraints from new data on peridotite massifs and xenoliths, 
Geochimica et Cosmochimica Acta, 70(17), 4528^550. 

Bitsch, B., A. Crida, A. S. Libert, and E. Lega (2013), Highly inclined and eccentric massive planets. 1. Planet-disc 
'\ntexd,ction^. Astronomy and Astrophysics, 555, 124. 

Blum, J., and G. Wurm (2008), The Growth Mechanisms of Macroscopic Bodies in Protoplanetary Disks, Annual 
Review of Astronomy and Astrophysics, 46(1), 21-56. 

Bottke, W. E, D. Nesvorny, D. Vokrouhlicky, and A. Morbidelli (2010), The Irregular Satellites: The Most Collision- 
ally Evolved Populations in the Solar System, The Astronomical Journal, 139(3), 994-1014. 

Bottke, W. E, D. Vokrouhlicky, S. Marchi, T. Swindle, E. R. D. Scott, and J. R. Weirich (2014), The Evolution of 
Giant Impact Ejecta and the Age of the Moon, 45th Lunar and Planetary Science Conference, 45, 1611. 

Brasser, R., K. J. Walsh, and D. Nesvorny (2013), Constraining the primordial orbits of the terrestrial planets. Monthly 
Notices of the Royal Astronomical Society, 433(A), 3417-3427. 

Briceno, C., A. K. Vivas, N. Calvet, L. Hartmann, R. Pacheco, D. Herrera, L. Romero, P. Berlind, G. Sanchez, J. A. 
Snyder, and P. Andrews (2001), The CIDA-QUEST Large-Scale Survey of Orion OBI: Evidence for Rapid Disk 
Dissipation in a Dispersed Stellar Population, Science, 291(5), 93-91. 

Bryden, G., M. Rozyczka, D. N. C. Lin, and P. Bodenheimer (2000), On the Interaction between Protoplanets and 
Protostellar Disks, The AstrophysicalJournal, 540(2), 1091-1101. 

Carballido, A., S. Eromang, and J. Papaloizou (2006), Mid-plane sedimentation of large solid bodies in turbulent 
protoplanetary discs. Monthly Notices of the Royal Astronomical Society, 373(4), 1633-1640. 

Carporzen, L., B. P. Weiss, L. T. Elkins-Tanton, D. L. Shuster, D. Ebel, and J. Gattacceca (2011), Erom the Cover: 
Magnetic evidence for a partially differentiated carbonaceous chondrite parent body, in Proceedings of the National 
Academy of Sciences, pp. 6386-6389. 

Castillo-Rogez, J. C., D. L. Matson, C. Sotin, T. V. Johnson, J. 1. Lunine, and P. C. Thomas (2007), lapetus’ geophysics: 
Rotation rate, shape, and equatorial ridge, ICARUS, 790(1), 179-202. 

Castillo-Rogez, J. C., T. V. Johnson, M. H. Lee, N. J. Turner, D. L. Matson, and J. Lunine (2009), 26A1 decay: Heat 
production and a revised age for lapetus, ICARUS, 204(2), 658-662. 

Chambers, J. E. (2001), Making More Terrestrial Planets, ICARUS, 152(2), 205-224. 

Chambers, J. E. (2007), On the stability of a planet between Mars and the asteroid belt: Implications for the Planet V 
hypothesis, ICARUS, 189(2), 386^00. 

Chambers, J. E. (2010), Planetesimal formation by turbulent concentration, ICARUS, 208(2), 505-517. 

Chambers, J. E. (2013), Late-stage planetary accretion including hit-and-run collisions and fragmentation, ICARUS, 
224(1), 43-56. 


19 



Chambers, J. E. (2014), Giant planet formation with pebble accretion, ICARUS, 233, 83-100. 

Chambers, J. E., and P. M. Cassen (2002), The effects of nebula surface density profile and giant-planet eccentricities 
on planetary accretion in the inner solar system, Meteoritics & Planetary Science, 37, 1523-1540. 

Chapman, C. R., B. A. Cohen, and D. H. Grinspoon (2007), What are the real constraints on the existence and 
magnitude of the late heavy bombardment?, ICARUS, 189(1), 233-245. 

Cohen, B. A., T. D. Swindle, and D. A. Kring (2000), Support for the Lunar Cataclysm Hypothesis from Lunar 
Meteorite Impact Melt Ages, Science, 290(5), 1754-1756. 

Crida, A., and A. Morbidelli (2007), Cavity opening by a giant planet in a protoplanetary disc and effects on planetary 
migration. Monthly Notices of the Royal Astronomical Society, 377(3), 1324-1336. 

Crida, A., Z. Sandor, and W. Kley (2008), Influence of an inner disc on the orbital evolution of massive planets 
migrating in resonance. Astronomy and Astrophysics, 483(1), 325-337. 

Cuk, M. (2012), Chronology and sources of lunar impact bombardment, ICARUS, 218(1), 69-79. 

Cuk, M., B. J. Gladman, and S. T. Stewart (2010), Constraints on the source of lunar cataclysm impactors, ICARUS, 
207(2), 590-594. 

Cuzzi, J. N., A. R. Dobrovolskis, and J. M. Champney (1993), Particle-gas dynamics in the midplane of a protoplane¬ 
tary nebula, ICARUS, 106, 102. 

Cuzzi, J. N., R. C. Hogan, and K. Shariff (2008), Toward Planetesimals: Dense Chondrule Clumps in the Protoplane¬ 
tary Nebula, The Astrophysical Journal, 687(2), 1432-1447. 

D’Angelo, G., and E. Marzari (2012), Outward Migration of Jupiter and Saturn in Evolved Gaseous Disks, The Astro- 
physicalJournal, 757(1), 50. 

Dauphas, N., and A. Pourmand (2011), Hf-W-Th evidence for rapid growth of Mars and its status as a planetary 
embryo. Nature, 473(1), 489^92. 

DeMeo, E. E., and B. Carry (2013), The taxonomic distribution of asteroids from multi-filter all-sky photometric 
surveys, ICARUS, 226(1), 723-741. 

DeMeo, E. E., and B. Carry (2014), Solar System evolution from compositional mapping of the asteroid belt. Nature, 
505(1), 629-634. 

Dominik, C., J. Blum, J. N. Cuzzi, and G. Wurm (2007), Growth of Dust as the Initial Step Toward Planet Eormation, 
Protostars and Planets V, pp. 783-800. 

Dubrulle, B., G. Morfill, and M. Sterzik (1995), The dust subdisk in the protoplanetary nebula, ICARUS, 114, 237-246. 

Elkins-Tanton, L. T., B. P. Weiss, and M. T. Zuber (2011), Chondrites as samples of differentiated planetesimals. Earth 
and Planetary Science Letters, 505(1), 1-10. 

Eernandez, J. A., and W. H. Ip (1984), Some dynamical aspects of the accretion of Uranus and Neptune - The exchange 
of orbital angular momentum with planetesimals, Icarus (ISSN 0019-1035), 58, 109-120. 

Eischer, R. A., and E. J. Ciesla (2014), Dynamics of the terrestrial planets from a large number of N-body simulations. 
Earth and Planetary Science Letters, 392, 28-38. 

Goldreich, P, and W. R. Ward (1973), The Eormation of Planetesimals, Astrophysical Journal, 183, 1051-1062. 

Gomes, R. S. (2003), The origin of the Kuiper Belt high-inclination population, ICARUS, 161(2), 404-418. 

Gomes, R. S., H. E Levison, K. Tsiganis, and A. Morbidelli (2005), Origin of the cataclysmic Late Heavy Bombard¬ 
ment period of the terrestrial planets. Nature, 435(1), 466^69. 


20 



Greenberg, R., J. F. Wacker, W. K. Hartmann, and C. R. Chapman (1978), Planetesimals to planets: Numerical 
simulation of collisional evolution, ICARUS, 55(1), 1-26. 

Guillot, T., S. Ida, and C. W. Ormel (2014), On the filtering and processing of dust by planetesimals 1. Derivation of 
collision probabilities for non-drifting planetesimals, arXiv.org, p. 7328. 

Haisch, K. E. J., E. A. Lada, and C. J. Lada (2001), Disk Erequencies and Lifetimes in Young Clusters, The Astrophys- 
ical Journal, 555(2), L153-L156. 

Halliday, A. N. (2008), A young Moon-forming giant impact at 70-110 million years accompanied by late-stage mix¬ 
ing, core formation and degassing of the Earth, Philosophical Transactions of the Royal Society A: Mathematical, 
Physical and Engineering Sciences, 566(1883), 4163^181. 

Hansen, B. M. S. (2009), Eormation of the Terrestrial Planets from a Narrow Annulus, The Astrophysical Journal, 
705(1), 1131-1140. 

Hopkins, P. E. (2014), Jumping the Gap: The Eormation Conditions and Mass Eunction of Pebble-Pile Planetesimals, 
arXiv.org, p. 2458. 

Ida, S., and D. N. C. Lin (2008), Toward a Deterministic Model of Planetary Eormation. IV. Effects of Type I Migration, 
The Astrophysical Journal, 675(1), 487-501. 

Ida, S., and J. Makino (1992), N-body simulation of gravitational interaction between planetesimals and a protoplanet. 
I - Velocity distribution of planetesimals, ICARUS, 96, 107-120. 

Ida, S., and J. Makino (1993), Scattering of planetesimals by a protoplanet - Slowing down of runaway growth, 
ICARUS, 106, 210. 

Ikoma, M., and H. Genda (2006), Constraints on the Mass of a Habitable Planet with Water of Nebular Origin, The 
Astrophysical Journal, 648{\), 696-706. 

Iwasaki, K., H. Emori, K. Nakazawa, and H. Tanaka (2002), Orbital Stability of a Protoplanet System under a Drag 
Eorce Proportional to the Random Velocity, Publications of the Astronomical Society of Japan, 54, 471^79. 

Izidoro, A., N. Haghighipour, O. C. Winter, and M. Tsuchida (2014), Terrestrial Planet Eormation in a Protoplanetary 
Disk with a Local Mass Depletion: A Successful Scenario for the Eormation of Mars, The Astrophysical Journal, 
7S2(1), 31. 

Jacobsen, S. B. (2005), The Hf-W Isotopic System and the Origin of the Earth and Moon, Annual Review of Earth and 
Planetary Sciences, 33, 531-570. 

Jacobson, S. A., and A. Morbidelli (2014), Lunar and terrestrial planet formation in the Grand Tack scenario, Phil 
Trans. R. Soc. A, 372, 0174. 

Jacobson, S. A., A. Morbidelli, S. N. Raymond, D. P. O’Brien, K. J. Walsh, and D. C. Rubie (2014), Highly siderophile 
elements in Earth’s mantle as a clock for the Moon-forming impact. Nature, 505(7494), 84-87. 

Jansson, K., and A. Johansen (2014), Eormation of pebble-pile planetesimals, arXiv.org, p. 2535. 

Jin, L., W. D. Arnett, N. Sui, and X. Wang (2008), An Interpretation of the Anomalously Low Mass of Mars, The 
Astrophysical Journal, 674(2), L105-L108. 

Johansen, A., and P. Lacerda (2010), Prograde rotation of protoplanets by accretion of pebbles in a gaseous environ¬ 
ment, Monthly Notices of the Royal Astronomical Society, 404, 475^85. 

Johansen, A., J. S. Oishi, M.-M. Mac Low, H. H. Klahr, T. Henning, and A. Youdin (2007), Rapid planetesimal 
formation in turbulent circumstellar disks. Nature, 448(1), 1022-1025. 

Johansen, A., J. Blum, H. Tanaka, C. Ormel, M. Bizzarro, and H. Rickman (2014), The multifaceted planetesimal 
formation process, arXiv.org, p. 1344. 


21 



Kita, N. T., G. R. Huss, S. Tachibana, Y. Amelin, L. E. Nyquist, and I. D. Hutcheon (2005), Constraints on the Origin 
of Chondrules and CAIs from Short-lived and Long-Lived Radionuclides, Chondrites and the Protoplanetary Disk, 
341, 558. 

Kley, W. (2000), The Orbital Evolution of Planets in Disks, Disks, 219, 69. 

Kley, W., and G. Dirksen (2006), Disk eccentricity and embedded planets. Astronomy and Astrophysics, 447(1), 369- 
377. 

Kobayashi, H., and N. Dauphas (2013), Small planetesimals in a massive disk formed Mars, ICARUS, 225(1), 122- 
130. 

Kobayashi, H., and H. Tanaka (2010), Lragmentation model dependence of collision cascades, ICARUS, 206(2), 735- 
746. 

Kokubo, E., and S. Ida (1995), Orbital evolution of protoplanets embedded in a swarm of planetesimals, ICARUS, 114, 
247-257. 

Kokubo, E., and S. Ida (1996), On Runaway Growth of Planetesimals, ICARUS, 123(1), 180-191. 

Kokubo, E., and S. Ida (1998), Oligarchic Growth of Protoplanets, ICARUS, 131(1), 171-178. 

Kretke, K. A., and H. E Levison (2014), Challenges in Eorming the Solar System’s Giant Planet Cores via Pebble 
Accretion, arXiv.org, p. 4430. 

Lambrechts, M., and A. Johansen (2012), Rapid growth of gas-giant cores by pebble accretion. Astronomy and Astro¬ 
physics, 544, 32. 

Lambrechts, M., and A. Johansen (2014), Eorming the cores of giant planets from the radial pebble flux in protoplan¬ 
etary discs. Astronomy and Astrophysics, 572, A107. 

Lambrechts, M., A. Johansen, and A. Morbidelli (2014), Separating gas-giant and ice-giant planets by halting pebble 
accretion. Astronomy and Astrophysics, 572, A35. 

Laskar, J. (1997), Large scale chaos and the spacing of the inner planets.. Astronomy and Astrophysics, 317, L75-L78. 

Lecuyer, C., P. Gillet, and E. Robert (1998), The hydrogen isotope composition of seawater and the global water cycle. 
Chemical Geology, 145(3), 249-261. 

Levison, H. E, and A. Morbidelli (2003), The formation of the Kuiper belt by the outward transport of bodies during 
Neptune’s migration. Nature, 426(6), 419^21. 

Levison, H. E, A. Morbidelli, C. Van Laerhoven, R. S. Gomes, and K. Tsiganis (2008), Origin of the structure of the 
Kuiper belt during a dynamical instability in the orbits of Uranus and Neptune, ICARUS, 196(1), 258-273. 

Levison, H. E, E. Thommes, and M. J. Duncan (2010), Modeling the Eormation of Giant Planet Cores. I. Evaluating 
Key Processes, The AstronomicalJournal, 139(A), 1297-1314. 

Levison, H. E, A. Morbidelli, K. Tsiganis, D. Nesvorny, and R. S. Gomes (2011), Late Orbital Instabilities in the 
Outer Planets Induced by Interaction with a Self-gravitating Planetesimal Disk, The Astronomical Journal, 142(5), 
152. 

Levison, H. E, M. J. Duncan, and D. M. Minton (2012), Modeling Terrestrial Planet Eormation with Lull Dynam¬ 
ics, Accretion, and Lragmentation, AAS/Division for Planetary Sciences Meeting Abstracts, 44. 

Levison, H. E, K. Kretke, and K. J. Walsh (2014), The Eormation of Terrestrial Planets from the Direct Accretion of 
Pebbles, AAS/Division for Planetary Sciences Meeting Abstracts, 46. 

Lin, D. N. C., and J. Papaloizou (1986), On the tidal interaction between protoplanets and the protoplanetary disk. Ill 
- Orbital migration of protoplanets. The Astrophysical Journal, 309, 846-857. 


22 



Lissauer, J. J. (1987), Timescales for planetary accretion and the structure of the protoplanetary disk, ICARUS, 69, 
249-265. 

Malhotra, R. (1995), The Origin of Pluto’s Orbit: Implications for the Solar System Beyond Neptune, Astronomical 
Journal v.llO, 110, 420. 

Mamajek, E. E. (2009), Initial Conditions of Planet Eormation: Lifetimes of Primordial Disks, in EXOPLANETS AND 
DISKS: THEIR FORMATION AND DIVERSITY: Proceedings of the International Conference. AIP Conference 
Proceedings, pp. 3-10, Department of Physics and Astronomy, University of Rochester, Rochester, NY 14627- 
0171. 

Marchi, S., W. E. Bottke, B. A. Cohen, K. Wunnemann, D. A. Kring, H. Y. McSween, M. C. de Sanctis, D. P. O’Brien, 
P. Schenk, C. A. Raymond, and C. T. Russell (2013), High-velocity collisions from the lunar cataclysm recorded in 
asteroidal meteorites. Nature Geoscience, 6(4), 303-307. 

Marchi, S., W. E. Bottke, L. T. Elkins-Tanton, M. Bierhaus, K. Wuennemann, A. Morbidelli, and D. A. Kring (2014), 
Widespread mixing and burial of Earth’s Hadean crust by asteroid impacts. Nature, 511(7), 578-582. 

Marty, B. (2012), The origins and concentrations of water, carbon, nitrogen and noble gases on Earth, Earth and 
Planetary Science Letters, 313, 56-66. 

Masset, E, and M. Snellgrove (2001), Reversing type II migration: resonance trapping of a lighter giant protoplanet. 
Monthly Notices of the Royal Astronomical Society, 320(4), L55-L59. 

Morbidelli, A., and A. Crida (2007), The dynamics of Jupiter and Saturn in the gaseous protoplanetary disk, ICARUS, 
797(1), 158-171. 

Morbidelli, A., J. E. Chambers, J. I. Lunine, J.-M. Petit, E Robert, G. B. Valsecchi, and K. E. Cyr (2000), Source 
regions and time scales for the delivery of water to Earth, Meteoritics & Planetary Science, 35, 1309-1320. 

Morbidelli, A., H. E. Levison, K. Tsiganis, and R. S. Gomes (2005), Chaotic capture of Jupiter’s Trojan asteroids in 
the early Solar System, Nature, 455(7041), 462^65. 

Morbidelli, A., K. Tsiganis, A. Crida, H. E. Levison, and R. S. Gomes (2007), Dynamics of the Giant Planets of the 
Solar System in the Gaseous Protoplanetary Disk and Their Relationship to the Current Orbital Architecture, The 
Astronomical Journal, 134(5), 1790-1798. 

Morbidelli, A., W. E. Bottke, D. Nesvorny, and H. E. Levison (2009), Asteroids were bom big, ICARUS, 204(2), 
558-573. 

Morbidelli, A., R. Brasser, R. S. Gomes, H. E. Levison, and K. Tsiganis (2010), Evidence from the Asteroid Belt for a 
Violent Past Evolution of Jupiter’s Orbit, The Astronomical Journal, 140(5), 1391-1401. 

Morbidelli, A., S. Marchi, W. E. Bottke, and D. A. Kring (2012a), A sawtooth-like timeline for the first billion years 
of lunar bombardment. Earth and Planetary Science Letters, 355, 144-151. 

Morbidelli, A., J. I. Lunine, D. P. O’Brien, S. N. Raymond, and K. J. Walsh (2012b), Building Terrestrial Planets, 
Annual Review of Earth and Planetary Sciences, 40, 251-275. 

Morbidelli, A., H. S. Gaspar, and D. Nesvorny (2014), Origin of the peculiar eccentricity distribution of the inner cold 
Kuiperbelt, ICARUS, 232, 81-87. 

Morishima, R., M. W. Schmidt, J. Stadel, and B. Moore (2008), Eormation and Accretion History of Terrestrial Planets 
from Runaway Growth through to Late Time: Implications for Orbital Eccentricity, The Astrophysical Journal, 
685(2), 1247-1261. 

Muralidharan, K., P. Deymier, M. Stimpfl, N. H. de Leeuw, and M. J. Drake (2008), Origin of water in the inner Solar 
System: A kinetic Monte Carlo study of water adsorption on forsterite, ICARUS, 198(2), 400^07. 

Murray-Clay, R., K. Kratter, and A. N. Youdin (2011), Core Accretion at Wide Separations: The Critical Role of Gas, 
AAS/Division for Extreme Solar Systems Abstracts, 2, 0804. 


23 



Nagasawa, M., D. N. C. Lin, and E. Thommes (2005), Dynamical Shake-up of Planetary Systems. 1. Embryo Trapping 
and Induced Collisions by the Sweeping Secular Resonance and Embryo-Disk Tidal Interaction, The Astrophysical 
Journal, 655(1), 578-598. 

Nakagawa, Y., C. Hayashi, and K. Nakazawa (1983), Accumulation of planetesimals in the solar nebula, ICARUS, 54, 
361-376. 

Nesvorny, D., D. Vokrouhlicky, and A. Morbidelli (2007), Capture of Irregular Satellites during Planetary Encounters, 
The AstronomicalJournal, 133(5), 1962-1976. 

Nesvorny, D., A. N. Youdin, and D. C. Richardson (2010), Eormation of Kuiper Belt Binaries by Gravitational Col¬ 
lapse, The Astronomical Journal, 140(3), 785-793. 

Nimmo, K, and T. Kleine (2007), How rapidly did Mars accrete? Uncertainties in the Hf W timing of core formation, 
ICARUS, 191(2), 497-504. 

O’Brien, D. P, A. Morbidelli, and H. E. Levison (2006), Terrestrial planet formation with strong dynamical friction, 
ICARUS, 184(\), 39-58. 

O’Brien, D. R, K. J. Walsh, A. Morbidelli, S. N. Raymond, and A. M. Mandell (2014), Water Delivery and Giant 
Impacts in the ’Grand Tack’ Scenario, ICARUS, 223, 74-84. 

Okuzumi, S., H. Tanaka, H. Kobayashi, and K. Wada (2012), Rapid Coagulation of Porous Dust Aggregates outside 
the Snow Line: A Pathway to Successful Icy Planetesimal Eormation, The Astrophysical Journal, 752(2), 106. 

Ormel, C. W., and H. H. Klahr (2010), The effect of gas drag on the growth of protoplanets. Analytical expressions 
for the accretion of small bodies in laminar disks. Astronomy and Astrophysics, 520, 43. 

Owen, T., and A. Bar-Nun (1995), Comets, Impacts, and Atmospheres, ICARUS, 116(2), 215-226. 

Papaloizou, J. C. B. (2003), Disc-Planet Interactions: Migration and Resonances in Extrasolar Planetary Systems, 
Celestial Mechanics and Dynamical Astronomy, 87(1), 53-83. 

Papaloizou, J. C. B., R. P. Nelson, and E. Masset (2001), Orbital eccentricity growth through disc-companion tidal 
interaction. Astronomy and Astrophysics, 366, 263-275. 

Pfalzner, S., M. Steinhausen, and K. Menten (2014), Short Dissipation Times of Proto-planetary Disks: An Artifact of 
Selection Effects?, The Astrophysical Journal Letters, 793(2), L34. 

Pierens, A., and R. P. Nelson (2008), Constraints on resonant-trapping for two planets embedded in a protoplanetary 
disc. Astronomy and Astrophysics, 482(1), 333-340. 

Pierens, A., and S. N. Raymond (2011), Two phase, inward-then-outward migration of Jupiter and Saturn in the 
gaseous solar nebula. Astronomy and Astrophysics, 533, 131. 

Pierens, A., S. Raymond, D. Nesvorny, and A. Morbidelli (2014), Outward migration of Jupiter and Saturn in 3:2 or 
2:1 resonance in radiative disks: implications for the Grand Tack and Nice models, arXiv.org, p. 543. 

Pollack, J. B., O. Hubickyj, P. Bodenheimer, J. J. Lissauer, M. Podolak, and Y. Greenzweig (1996), Eormation of the 
Giant Planets by Concurrent Accretion of Solids and Gas, ICARUS, 124(1), 62-85. 

Raymond, S. N., and A. Morbidelli (2014), The Grand Tack model: a critical review, eprint arXiv:1409.6340. 

Raymond, S. N., T. R. Quinn, and J. 1. Lunine (2004), Making other earths: dynamical simulations of terrestrial planet 
formation and water delivery, ICARUS, 168(1), 1-17. 

Raymond, S. N., D. P. O’Brien, A. Morbidelli, and N. A. Kaib (2009), Building the terrestrial planets: Constrained 
accretion in the inner Solar System, ICARUS, 203(2), 644-662. 

Raymond, S. N., E. Kokubo, A. Morbidelli, R. Morishima, and K. J. Walsh (2013), Terrestrial Planet Eormation at 
Home and Abroad, arXiv.org, p. 1689. 


24 



Ryder, G. (1990), Lunar samples, lunar accretion and the early bombardment of the moon, EOS (ISSN 0096-3941), 
77,313. 

Ryder, G. (2002), Mass flux in the ancient Earth-Moon system and benign implications for the origin of life on Earth, 
Journal of Geophysical Research (Planets), 707(E), 5022. 

Safronov, V. S. (1972), Evolution of the protoplanetary cloud and formation of the earth and planets.. Evolution of the 
protoplanetary cloud and formation of the earth and planets. 

Safronov, V. S., and E. V. Zviagina (1969), Relative sizes of the largest bodies during the accumulation of planets, 
ICARUS, 70(1), 109-115. 

Snellgrove, M. D., J. C. B. Papaloizou, and R. R Nelson (2001), On disc driven inward migration of resonantly coupled 
planets with application to the system around GK76, Astronomy and Astrophysics, 374, 1092-1099. 

Strom, R. G., R. Malhotra, T. Ito, E. Yoshida, and D. A. Kring (2005), The Origin of Planetary Impactors in the Inner 
Solar System, Science, 309(5), 1847-1850. 

Taylor, D. J., K. D. McKeegan, and T. M. Harrison (2009), Lu-Hf zircon evidence for rapid lunar differentiation. Earth 
and Planetary Science Letters, 279(3), 157-164. 

Tera, K, D. A. Papanastassiou, and G. J. Wasserburg (1974), Isotopic evidence for a terminal lunar cataclysm. Earth 
and Planetary Science Letters, 22(1), 1-21. 

Thommes, E., M. Nagasawa, and D. N. C. Lin (2008), Dynamical Shake-up of Planetary Systems. 11. N-Body Sim¬ 
ulations of Solar System Terrestrial Planet Eormation Induced by Secular Resonance Sweeping, The Astrophysical 
Journal, 676(1), 728-739. 

Touboul, M., T. Kleine, B. Bourdon, H. Palme, and R. Wieler (2007), Late formation and prolonged differentiation of 
the Moon inferred from W isotopes in lunar metals. Nature, 450(1), 1206-1209. 

Tsiganis, K., R. S. Gomes, A. Morbidelli, and H. E. Levison (2005), Origin of the orbital architecture of the giant 
planets of the Solar System, Nature, 435(1), 459^61. 

Villeneuve, J., M. Chaussidon, and G. Libourel (2009), Homogeneous Distribution of 26A1 in the Solar System from 
the Mg Isotopic Composition of Chondrules, Science, 325(5), 985. 

Wada, K., H. Tanaka, T. Suyama, H. Kimura, and T. Yamamoto (2008), Numerical Simulation of Dust Aggregate Col¬ 
lisions. II. Compression and Disruption of Three-Dimensional Aggregates in Head-on Collisions, The Astrophysical 
Journal, 677(2), 1296-1308. 

Wada, K., H. Tanaka, T. Suyama, H. Kimura, and T. Yamamoto (2009), Collisional Growth Conditions for Dust 
Aggregates, The Astrophysical Journal, 702(2), 1490-1501. 

Wada, K., H. Tanaka, S. Okuzumi, H. Kobayashi, T. Suyama, H. Kimura, and T. Yamamoto (2013), Growth efficiency 
of dust aggregates through collisions with high mass ratios. Astronomy and Astrophysics, 559, A62. 

Walker, R. J. (2009), Highly siderophile elements in the Earth, Moon and Mars: Update and implications for planetary 
accretion and differentiation, Chemie der Erde - Geochemistry, 69, 101-125. 

Walsh, K. J., and H. E. Levison (2014), Eormation of the Terrestrial Planets from an Annulus, AAS/Division of Dy¬ 
namical Astronomy Meeting, 45. 

Walsh, K. J., and A. Morbidelli (2011), The effect of an early planetesimal-driven migration of the giant planets on 
terrestrial planet formation. Astronomy and Astrophysics, 526, 126. 

Walsh, K. J., A. Morbidelli, S. N. Raymond, D. P. O’Brien, and A. M. Mandell (2011), A low mass for Mars from 
Jupiter’s early gas-driven migration. Nature, 475(1), 206-209. 


25 



Walsh, K. J., A. Morbidelli, S. N. Raymond, D. R O’Brien, and A. M. Mandell (2012), Populating the asteroid belt 
from two parent source regions due to the migration of giant planets—"The Grand Tack", Meteoritics & Planetary 
Science, 47(1), 1941-1947. 

Ward, W. R. (1980), Scanning Secular Resonances: a Cosmogonical Broom?, LUNAR AND PLANETARY SCIENCE 
XL 11, 1199-1201. 

Weidenschilling, S. J. (1977), Aerodynamics of solid bodies in the solar nebula. Monthly Notices of the Royal Astro¬ 
nomical Society, 180, 57-70. 

Weidenschilling, S. J. (1980), Dust to planetesimals - Settling and coagulation in the solar nebula, ICARUS, 44, 172- 
189. 

Weidenschilling, S. J. (2011), Initial sizes of planetesimals and accretion of the asteroids, ICARUS, 214(2), 671-684. 

Weidenschilling, S. J., and J. N. Cuzzi (1993), Formation of planetesimals in the solar nebula. Protostars and planets 
III, pp. 1031-1060. 

Weiss, B. R, J. Gattacceca, S. Stanley, P. Rochette, and U. R. Christensen (2010), Paleomagnetic Records of Meteorites 
and Early Planetesimal Differentiation, Space Science Reviews, 752(1), 341-390. 

Wetherill, G. W. (1985), Occurrence of giant impacts during the growth of the terrestrial planets. Science (ISSN 0036- 
8075), 228, 877-879. 

Wetherill, G. W. (1991), Why Isn’t Mars as Big as Earth?, Abstracts of the Lunar and Planetary Science Conference, 
22, 1495. 

Wetherill, G. W., and G. R. Stewart (1989), Accumulation of a swarm of small planetesimals, ICARUS, 77, 330-357. 

Wetherill, G. W., and G. R. Stewart (1993), Eormation of planetary embryos - Effects of fragmentation, low relative 
velocity, and independent variation of eccentricity and inclination, ICARUS, 106, 190. 

Whipple, E. L. (1972), On certain aerodynamic processes for asteroids and comets. From Plasma to Planet, p. 211. 

Yin, Q., S. B. Jacobsen, K. Yamashita, J. Blichert-Toft, P. Telouk, and E. Albarede (2002), A short timescale for 
terrestrial planet formation from Hf-W chronometry of meteorites. Nature, 475(6901), 949-952. 

Youdin, A. N. (2010), Erom Grains to Planetesimals, EAS Publications Series, 41, 187-207. 

Youdin, A. N., and E. H. Shu (2002), Planetesimal Eormation by Gravitational Instability, The Astrophysical Journal, 
580(1), 494-505. 

Zhou, J.-L., D. N. C. Lin, and Y.-S. Sun (2007), Post-oligarchic Evolution of Protoplanetary Embryos and the Stability 
of Planetary Systems, The AstrophysicalJournal, 666(1), 423^35. 

Zsom, A., C. W. Ormel, C. Guttler, J. Blum, and C. P Dullemond (2010), The outcome of protoplanetary dust growth: 
pebbles, boulders, or planetesimals? II. Introducing the bouncing barrier. Astronomy and Astrophysics, 513, 57. 


26 



